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We present a numerical scheme to study the dynamics of slow light and light storage in an 
eletromagnetically induced transparency (EIT) medium at finite temperatures. Allowing for the 
moitonal coupling, we derive a set of coupled Schrodinger equations describing a boosted closed 
3-level EIT system according to the principle of Galilean relativity The dynamics of a uniformly 
moving EIT medium can thus be determined by numerically integrating the coupled Schrodinger 
equations for atoms plus one ancillary Maxwell- Schrodinger equation for the probe pulse. The 
central idea of this work rests on the assumption that the loss of ground-state coherence at finite 
temperatures can be ascribed to the incoherent superposition of density matrices representing the 
EIT systems with various velocities. Glose agreements are demonstrated in comparing the numerical 
results with the experimental data for both slow light and light storage. In particular, the distinct 
characters featuring the decay of ground-state coherence can be well verified for slow light and light 
storage. This warrants that the current scheme can be applied to determine the decaying profile of 
the ground-state coherence as well as the temperature of the EIT medium. 

PACS numbers: 42.50.Gy, 32.80.Qk, 02.60.Jh 



I. INTRODUCTION 



The effect of electromagnetically induced transparency 
(EIT) is a nonlinear optical phenomenon which renders 
an opaque medium transparent at a certain frequency by 
exciting it with an electromagnetic field [l] . This effect 
can occur generally in a 3-level atomic system, where 
the atomic states are coherently prepared by external 
laser fields. For example, in a A- type system (see Fig.l), 
there are two dipole- allowed transitions |1) ^ |3) and 
1 2) ^ 1 3) which are excited by a weak probing field, 
and a strong coupling field, respectively. When EIT oc- 
curs, destructive interference among different transition 
pathways suppresses the transition probability between 
|1) ^ 1 3), leading to a spectrally sharp dip in the ab- 
sorption spectrum. The corresponding steep dispersion 
within the transparency window results in a large reduc- 
tion in the group velocity of light. As the steepness de- 
pends on the intensity of coupling field, EIT thus provides 
an effective and convenient mechanism for slowing down 
the light in a controllable fashion. The slow light (SL) 
arising from the EIT effect greatly enhances the nonlin- 
ear susceptibility and makes the low-light level nonlinear 
optics possible [2 . The first experimental demonstration 
of SL produced by EIT was made with high-power pulsed 
lasers interacting with a Sr vapor p|||4) . In 1999, a dra- 
matic reduction of the group velocity down to 17m/s was 
demonstrated by Hau et at . by using a Bose-Einstein 
condensate of Na atoms [5 . An ultraslow group velocity 
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of 8m/s were later observed in a buffer-gas cell of hot Rb 
atoms by Budke et al . [6 . 

In a lossless, passive sample, the reduction of the prop- 
agation velocity implies a temporary transfer of electro- 
magnetic excitations into the medium. With EIT, the 
light pulse can even be completely stopped and stored 
in the medium by adiabatically switching off the cou- 
pling field and subsequently retrieved from the medium 
by the reverse process while the probe pulse is entirely 
within the sample making its way through [11 [6l [71 [11] 
. Such storage and retrieval of photonic information are 
essentially reversible as they are resulted by the coher- 
ent transfer of the quantum state of light into the quan- 
tum coherence of the two ground states |1) and |2). The 
light storage (LS) was experimentally demonstrated by 
several groups with various schemes [8] fTTHT3] , which 
promises to be applied in the processing of quantum in- 
formation, especially in the implementation of quantum 
storage devices, logical gates and generation of photonic 
qubits [HHTS] . 

So far, except only a few experiments were demon- 
strated in atomic Bose-Einstein condensates at nearly 
zero temperature ^ (8j , most experimental studies for 
EIT were carried out at finite temperatures, using ei- 
ther hot atoms at room temperatures |9] or laser- 
cooled atoms at temperatures about a few hundreds of 
microKelvin [T] . It is well-known that at finite tempera- 
tures, the Doppler-broadened medium inevitably imposes 
a serious limit since the Doppler shifts caused by the 
atomic motion introduce a randomization in the effective 
laser detunings over the ensemble of atoms in the sample 
[To] . Experimental evidences indicate that even in the 
laser-cooled atoms, the decoherence due to the atomic 
motion still can not be ignored readily, especially when 
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the applied external fields are in a counter-propagating 
geometry [iTl [18] . Phenomenologically, the effect of 
Doppler broadening can be addressed by including a re- 
laxation term in the motion equation of pi2. Here pi2 
denotes the off-diagonal element of the density operator 
which describes the quantum coherence of the two ground 
states |1) and |2). Accordingly, a decay constant 7 is thus 
introduced to account for the relaxation of ground-state 
coherence. It should be noted that 7 is by no means 
a universal constant which depends on the temperature 
as well as on the experimental parameters such as the 
amplitudes of the driving fields and the particle density 
of the medium. Normally, 7 can only be determined by 
numerically fitting upon the experimental data. Thus 
one has to recompute the numerical value of 7 once the 
experimental parameters are altered in the new runs of 
experiment. In order to determine the dynamical proper- 
ties of the system without introducing 7, it is desirable to 
develop a numerical scheme which is tractable and acces- 
sible to both theorists and experimentalists, and in this 
paper we address such an issue. 

In general, the atomic dynamics of the EIT medium 
can be approached by employing the formalism of master 
equation in the Linblad form, 

'^ = i[H,m]+vim) (1) 

where p {t) is the density matrix of the 3-level atom, H 
is the total Hamiltonian, and V {p (t)) is the dissipator of 
the master equation through which the relevant decaying 
processes such as the spontaneous emission and absorp- 
tion can be specified. Essentially, the medium is envis- 
aged as a single, huge and stationary, 3-level "atom", so 
that only the internal atomic degrees of freedom are con- 
sidered in the last equation. To include the effects of non- 
relativistic thermal random motion of atoms, a plausible 
way is to add a Doppler energy term fik • v directly to the 
Hamiltonian in eq.Q, where k is the wave vector of the 
applied laser field, and v is the relative velocity between 
the atom and the light source. On account of the ther- 
mal motion of atoms, v is randomized, and therefore the 
final solution of eq.Q has to be determined statistically. 
Denoting p*^^^ as the solution of eq.([I]) with a given v, 
then the desired density matrix can be obtained by tak- 
ing average over all possible velocities which are typically 
described by Maxwell-Boltzmann distribution when the 
gas is at thermal equilibrium. It is expected that such an 
averaging process smears the atomic coherence and thus 
contributes to the overall effect of decoherence [10] . 

We note that, for taking the effect of Doppler shift into 
account the above formula is valid only when the "atom" 
is stationary or, equivalently, when the frame of reference 
is fixed on the "atom" . As we aim to obtain the statisti- 
cally averaged density matrix over an ensemble of identi- 
cal EIT systems moving with different velocities, a partic- 
ular frame of reference must be specified in which each 
p^^^ can be solved and properly weighted on the com- 
mon ground. This would involve the transformation of 



motion equations under the principle of Galilean relativ- 
ity. However, in the presence of driving fields, the atoms 
are by no means in constant motion as they do exchange 
momenta with the light fields by absorbing or emitting 
photons from time to time. Furthermore, the atoms are 
liable to couple to vacuum via spontaneous emission of 
photons to all directions, giving rise to the random re- 
coil of the atoms. These circumstances suggests that the 
moving atoms do not serve as an ideal frame of reference 
in our scheme. Alternatively, we may seek to solve the 
matrix elements of all p^^^ {t) and carry out the ensem- 
ble average in the laboratory frame [19] . In doing so, 
we might have to consider the sample g as as a lump of 
medium rather than a single "atom". Accordingly, it is 
more convenient to work in the framework of Schrodinger 
equation instead of the master equation. There are two 
reasons: the gauge invariance of Schrodinger equation 
under the Galilean transformation are well understood. 
Second, the kinetic energy can be naturally included in 
the Schrodinger equation, so that the motional coupling, 
namely, the coupling between the external degrees of free- 
dom and the internal degrees of freedom, can be restored 
and properly dealt with, although they are normally over- 
whelmed by the light-atom interaction. The details of 
the derivation of the motion equations will be described 
later. 

The organization of this paper is as follows. In Sec. II, 
we derive the general forms of the motion equations of 
the atomic medium under arbitrary Galilean transforma- 
tion boosted along the z-direction. The numerical results 
are presented in Sec. HI. Comparison with experimen- 
tal data are made. Finally, some concluding remarks are 
given in Sec. IV. 

II. FORMALISM 

We consider a medium consisting of A-type 3-level 
atoms with two met ast able ground states as shown in 
Fig. 1. The atoms in this medium are all noninteract- 
ing and excited by two laser fields. The probe pulse, 
which drives the |1) ^ |3) transition is characterized by 
a central frequency Up and a wave vector kp. On the 
other hand, the |2) ^ |3) transition is driven by another 
laser field with frequency ujc and wave vector kc. For 
all practical purposes, the probe and couple fields can be 
applied with different relative orientation, depending on 
the experimental applications. For simplicity, we shall 
assume in the following derivations that the probe field 
propagates along the z direction and the coupling field 
propagates with an angle 6 with respect to the z-axis 
and the 3-level atoms are in a cigar-shaped trap which 
can be considered as an one dimensional system. Now 
the Hamiltonian of the system can be given by 

H = Ho + Hi, (2) 

where 
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and consequently, we obtain 
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dt 



ih 



d^2 
dt 



2m dz^ 



2m dz'^ 



(8) 



^2 + -!^:e-<^=^-'^»*)^3, 



(9) 



FIG. 1: (Color online) The energy levels of a 3- level atom of 
A- type. The transitions |1) ^ |3) and |2) ^ |3) are driven by 
laser fields with Rabi frequencies Vtp and Qc correspondingly. 
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is the unperturbed part and 
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describe the interaction between EM field and the atom 
in the dipole approximation and rotating-wave approxi- 
mation. Here, Ej is the energy of the electronic level 
kc = kc cos 6 and and ftp denote the Rabi frequencies 
for the transitions |2) ^ |3) and |1) ^ |3) respectively. 
In the current problem, we assume that Qp = ^p{z^ t) 
and Qc = (^)- 

In the continuum limit, a 3- level atomic medium can 
be described by a 3-component spinor field. 
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where ?/^j(r,t) is the atomic field operator which annihi- 
lates an atom in the internal state \j) that is positioned 
at z. In terms of the atomic field operators, the energy 
functional of the above EIT system is given by 
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The motion equations for all tpj (r, t) can be derived from 
the Hartree variational principle, namely, 



(i = 1,2,3) 



(7) 



Here we have ignored the transverse motion of the atoms 
in the x?/-plane since the Gaussian probe pulse propa- 
gates along the z-direction only. 

It is convenient to describe the atomic properties by 
means of the local density operator p whose matrix el- 
ements are defined as bilinear products of atomic fields, 
i.e., pij = il^iil^"^. Accordingly, the z-th diagonal matrix 
element corresponds to the density of atoms in the state 
In the absence of any decaying process, the total den- 
sity may be normalized to / dz = NL^ where N 
is total particle number and L is the length of the sys- 
tem. The quantum coherence between the states \i) and 
\j) is represented by the off-diagonal matrix element pij. 
In particular, the matrix element, pi2, is of central im- 
portance in the current problem, which dominates the 
dynamics of the storage and retrieval of the probe light 
pulse. Additionally, the coherence between |1) and |3) 
determines the propagation of the probe field inside the 
medium. In the slowly varying envelope approximation 
[20] , the dynamics of the probe field is governed by the 
Maxwell- Schrodinger equation 
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where r] = SXj^MaT /An with A/" being the number per 
unit length of the medium), T the spontaneous decay 
rate of |3), a the branch ratio of the decay from |3) to 
and Xl being the wavelength of the laser field. 
Note that eqs.(l8])-( 10) can be interpreted as the gov- 



erning equations of a moving continuous medium seen 
by a co-moving observer, although they take the form 
of coupled Schrodinger equations. In the nonrelativis- 
tic quantum mechanics, the Galilean relativity ensures 
that the Schrodinger equation is form- invariant under the 
Galilean transformation, r ^ r — vt, t ^ for any two 
frames of reference that are in relative uniform transla- 
tional motion with a velocity v [21 . To be specific, let 
us assume that the boosted system moves with a veloc- 
ity V = vz relative to the laboratory frame. Accordingly, 
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the Galilean principle of relativity demands the gauge 
dependence of two different but equivalent quantum me- 
chanical states 

* (z, t) = (z^ t') e-<-^'V2-m^.)/^^ ^^2) 

where (z, t) and {z' ^ t') denote the space-time coordinates 
in the laboratory frame and the boosted frame respec- 
tively, which are transformed hy z' = z — vt^ t' = t. 
Moreover, we require that the pulse profiles in different 
frames are related by, 

(z, t) = % (z^ t') = %{z- vt, t) , (13) 

and thus eqs.([8|-([To]) become 
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+ ^ne{t)e<^^^--^')4''\z-vt,t), 

where tpj^^ {j = 1,2,3) denote the atomic fields in the 
medium boosted with velocity v and Uj = Ej/h. 

To solve eqs.([l4|-([l6| in a more efficient manner (less 
grid points and higher accuracy), we extract a fast oscil- 
lating phase (both spatially and temporarily) from each 
by writing 



spontaneous decay of the excited level |3) into account, 

(v) 

we then obtain the following motion equations for ^ in 
the laboratory frame. 
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where Ap = ujs — oJi — ojp^ Ac = ujs — uj2 — are the 
detunings of probe and coupling lasers respectively, and 

= (^kp — kc^ V + h{kp — kcY /2m denotes the motion- 
induced frequency shift. In eq.([20|, the spontaneous de- 
cay of the excited level |3) occurring with a rate T is in- 
cluded phenomenologically. Note that the term, kpV and 
kcV entering the right-hand-side of eq.(19) and eq.(20), 
are the nonrelativistic Doppler shift which turn out to 
be the detunings for the probe and coupling laser respec- 
tively in the laboratory frame. In addition, it is note- 
worthy to point out that h{kp — kcY /2m represents the 
frequency shift caused by the recoil of atom. 

The dephasing owing to the thermal motion of atoms 
is incorporated in the system dynamics by replacing pai 
with (p3i)^ in the right-hand-side of eq.(ll), namely. 
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where ^^-^^ represents the slowly varying part for 
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Substituting eq.(17) into eqs.( 14)-( 16) and taking the 
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2kBTi^ 



-mv^/2kBT 



the thermally averaged atomic coherence between |1) and 
1 3), which is taken over the Maxwell-Boltzmann distribu- 
tion of velocity at a given temperature T. 
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Finally, we give a brief account for the numerical 
method which we have employed to integrate eqs.(18)- 
(21). Since we have used very high resolution to re- 



solve the fine structures of the wave functions and the 
light pulse during their time evolutions, the commonly 
used second-order Crank-Nicolson method turns out to 
be inefficient in the current problem despite that it is un- 
conditionally stable. Owing to its implicitness in time, 
the Crank-Nicolson method would require an exceed- 
ingly large resultant matrix for the four coupled equa- 
tions, eqs.(18)-(21 ), when a high spatial resolution is de- 



manded. With this regard, explicit methods are more ef- 
ficient for the current problem. Here we use the method 
of lines with spatial discretization by highly accurate 
Fourier pseudospectral method and time integration by 
adaptive Runge-Kutta method of orders 2 and 3 (RK23), 
such that the accuracy in time is at least of second order, 
and the accuracy in space is of exponential order. In the 
following numerical computations, we will first determine 
(j)^^^ (z, t) by numerical integration. Having obtained all 
(f)^^^ (z, t), we then calculate (pai (z, t))^ to determine the 



profile of Vtp at the instant t. 



III. RESULTS AND DISCUSSIONS 

In this section, we apply the aforementioned scheme to 
simulate the effect of Doppler broadening for both SL and 
LS, and compare the numerical results with the experi- 
mental data. Before we present the results, some remarks 
are given. To begin with, let us consider the limiting 
case with kp = kc and 6> = 0, where the shift in the 



right-hand-side of eq.(19) can be exactly cancelled out. 



indicating that the EIT dynamics can be robust against 
the effect of Doppler broadening when probe and cou- 
ple are co-propagating. However, we note that in most 
experiments, kp differs from kc in the transition scheme, 
or a small is required to separate the weak probe and 
strong coupling fields in the detection scheme. This sug- 
gests that the effects of Doppler broadening is inevitable 
in practical applications, despite that they are essentially 
small in the co-propagating situation. On the contrary, 
setting kp = kc and = it, Ay is maximized for all v and 
can not be eliminated in any case and thus a substantial 
loss of ground-state coherence is expected in the counter- 
propagating geometry. Regarding the fact that the deco- 
herence is much pronounced in the latter case, plus that 
the counter-propagating EIT is the central mechanism 
for generating stationary light pulses [17^ , here we re- 
strict our attention on the counter-propagating EIT. As 
a matter of fact, in our co-propagating EIT experiment 
in which the effect of the thermal motion of the cold 
atoms is minimized, the decoherence rate is found to be 
less than 10"^r (F = 27r x 5.9 MHz), indicating that the 
relaxations caused by some intrinsic effects, such as the 
stray magnetic field, collisions, and linewidth and fre- 
quency fluctuation of the laser fields, etc., are very small 



which add up to give a decoherence rate less than 10~^F. 
Such a small decoherence rate suggests that the decay of 
the output probe fields in our counter-propagating ex- 
periment are mainly from the thermal motion. In fact, 
adding such small decoherence rate in our calculation 
changes the numerical outcome very little and therefore 
we shall not consider the intrinsic decoherence in the nu- 
merical simulations. 

The comparisons between experimental data and nu- 
merical results of SL are made in Fig. 2 for three dif- 
ferent experiments. We carried out the measurements 
in a cigar-shaped cloud of laser-cooled Rb atoms [22] in 
which the probe pulse and coupling field are counter- 
propagating along the major axis of the cloud. The 
Rabi frequency Qc and the optical density OD (= 2Lr]/r) 
for the three samples are estimated to be Qc = O.SSF, 
OD = 30, Qc = 0.69F, OD = 41 and Qc = 0.77F, 
OD = 48 respectively. Both Clc and OD are estimated 
by the method described in fT7| and have an uncertainty 
of about ±5%. The probe pulses used in the experiment 
are sufficiently weak such that the corresponding Rabi 
frequency in the calculation does not affect the prediction 
of the output probe pulse. In the numerical simulations 
with 128 velocity groups, we set Qc = 0.825F, OD = 30, 
T = 290 /iK, Qc = 0.665F, OD = 41.3, T = 305 fiK and 
Qc = 0.75r, OD = 48, T = 280 /iK to get the best fit for 
the three experiments of SL shown in Fig.2(a)-2(c), re- 
spectively. It should be mentioned that the temperature 
of the SL experiment in Fig. 2 (a) has been determined by 
another different numerical method to give the value of 
T = 296 fiK [23] , which is very close to our prediction. 
The experimental data and numerical results of storage 
and retrieval of a light pulse are shown in Fig. 3. In the 
simulations of Fig. 3, we only adjust the temperature to 
T = 240 /iK while keep the values of Qc and OD in 
Fig. 2 (a) unchanged to get the best fit. The discrepancies 
in the temperatures obtained from the numerical simu- 
lations in Fig. 2 (a) and Fig. 3 are very reasonable as com- 
pared with the expected temperature and fluctuations of 
the laser-cooled Rb atoms, yet Clc and OD so obtained 
are in good agreement with the experimental parameters. 
Therefore, we have quantitatively demonstrated the va- 
lidity and accuracy of the numerical method presented 
in this work. In addition, we find this method can be 
used to determine the temperature along the major axis 
of the cigar-shaped atom cloud which we were not able 
to measure previously. 

The influence of the thermal effect is readily revealed 
by the considerably diminished intensity of the output 
probe beam. In Fig. 4(a), the intensity profiles of the 
output probe pulse at various temperatures are shown 
for SL. As expected, the output intensity decreases when 
the temperature increases. A simple explanation for this 
is that the higher the temperature, the broader the veloc- 
ity distribution is. Thus the effective two-photon detun- 
ing (the averaged Doppler shift) becomes larger and the 
output probe pulse is significantly suppressed. Fig. 4(b) 
shows that the peak of the output probe decays exponen- 
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tially as a function of temperature. 
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FIG. 2: (Color online) The best fitting of the numerical sim- 
ulation with (a) Qc = 0.825r, OD = 30, T = 290 /xK; (b) 
= 0.665r, OD = 41.3, T = 305 /xK; (c) Qc = 0.75r, 
OD = 48, T = 280 /iK for SL. The squares and circles rep- 
resent the intensities of the experimentally measured input 
probe pulse and output probe pulse, respectively. The input 
probe pulse is scaled down by a factor of 0.2 in (a), 0.01 in 
(b) and 0.04 in (c), respectively. The dashed and solid lines 
represent the nurmerical simulations for the intensities of the 
input probe pulse and the output probe pulse, respectively. 



For all practical purposes, it is more instructive to ex- 
amine the ground-state coherence, {p2i)t ^ rather than 
the output intensity of the probe beam, since there are 
no light fields during the process of storage, and once the 
stored pulse is retrieved, it restores SL again. Although 
(P2i)t ^ot be measured directly, its magnitude deter- 
mines the intensity of the retrieved probe pulse when the 
coupling field is turned off. Very recently, in analyzing 
the feasibility of measuring the ground-state coherence 
in an EIT, Zhao et al. [19] suggested that (p2i)t decays 
like a Gaussian function during LS. On this basis, it can 
be further shown that 



{p2i [z, t))j. = {p2i{z, t = 0))^ e 



-{kvs) J2ikz 



(22) 



where Vs = ^j2kBT /m is the one-dimensional root-mean- 
square velocity and k = c~^{ujs — ui). Now let us con- 
sider the case of the storage time of 1.4 jas for various 
temperatures T = 100, 300, 500 jaK. In the calculation. 
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FIG. 3: (Color online) The best fitting of the numerical sim- 
ulation with Qc = 0.825r, OD = 30, T = 240 /xK for LS. 
The squares and circles are the experimentally measured in- 
tensities of the input probe pulse and output probe pulse, 
respectively. The dashed line and solid line represent the the 
nurmerical simulations for the intensities of the input probe 
pulse and output probe pulse, respectively. The input probe 
pulse is scaled down by a factor of 0.05. 



the probe pulse enters the medium under Qc = 0.75r 
and OD = 32 both of which do not affect the decay be- 
havior during the storage; the input probe pulse and the 
timing of switching off the coupling field are the same as 
those shown in Fig. 3. For such a storage time and atom 
temperatures the atomic thermal motion is expected to 
well smear the quantum memory of the probe pulse which 
is stored as the ground-state coherence. Because eq.(22) 
is a function of z and t and the time-dependence only 
comes from the Gaussian function, to eliminate the z 
dependence of the ground-state coherence, we plot the 



function R21 {t^T) = |J (p2i {z,t))j.dz\ in Fig. 5. We see 
that the R21 (t, T) indeed decays Gaussian-like with time 
as predicted in \T9^ . Now let us denote r as the 1/e 
width of the Gaussian function. Accordingly, r of the 
fitted Gaussian function of the ground-state coherence 
in the simulation is found to be close to the value pre- 
dicted by eq.(22). Furthermore, we have verified that 
r oc T~^/^, and the close agreement with the theoreti- 
cal predictions of eq.(22) are also shown in Fig. 5. Since 



the thermal motion would randomize the spatial profile 
of the ground-state coherence during the storage, the in- 
tensity of the retrieved probe pulse is thus much smaller 
than the stored one as shown in Fig. 3. 

The ground-state coherence for SL can be studied in 
a similar manner. It is expected that the ground-state 
coherence decays more slowly in SL than in LS, since the 
continuous optical excitations from the coupling field can 
retard the loss of atomic coherence led by the randomiza- 
tion of atoms' thermal motion. We simulate the process 
of SL with ftc = 0.75r, and OD = 400 at various temper- 
atures. Here we have chosen a very large optical density 
OD = 400 to ensure that the probe pulse can stay in the 
medium for a sufficiently long period. We plot R21 (t, T) 



7 



as a function of time for various temperatures after the 
probe pulse has entirely merged into the medium, and 
the numerical results suggest that R21 (t, T) decays ex- 
ponentially with a rate n. In the following, we apply some 
previous theoretical results based on the steady-state so- 
lution of the optical-Bloch equation to derive an analytic 
estimate of n. For simplicity, we shall assume that the 
EIT transparency bandwidth is much larger than the fre- 
quency bandwidth of the probe pulse, such that the de- 
cay of the probe pulse or the ground-state coherence is 
mainly determined by the absorption in the center of the 
transparency window of the EIT spectrum. Because all 
decoherence mechanisms other than the atomic thermal 
motion are neglected, the steady-state solution of psi of 
a rest EIT medium is given by [24] 



(a) 



O 





- inputxO.15 

-- -output T=1 00 |iK 
-- -output T=1 50 |xK 

- - outputT=200|iK 

- - output T=250|iK 
output T=300|iK 







Time [ |is] 



O output peak 
— fitting curve 




200 

Temperature [ |iK] 



FIG. 4: (Color online) The best fitting of the numerical sim- 
ulation with = 0.825r, OD = 30, T = 240 /xK for LS. 
The squares and circles are the experimentally measured in- 
tensities of the input probe pulse and output probe pulse, 
respectively. The dashed line and solid line represent the the 
nurmerical simulations for the intensities of the input probe 
pulse and output probe pulse, respectively. The input probe 
pulse is scaled down by a factor of 0.05. 



Psi 



A, - A„ 



n2/2 + 2[Aj,-ir/2][Ae- Ap]' 



(23) 



Given that the probe and coupling fields are both reso- 
nant with the atomic transition frequencies in the labo- 
ratory frame, thus for an atom moving with a velocity 



In the presence of a strong coupling field, Qc ^ kvg, the 
imaginary part of eq.(23) can be approximated by 



Im 



(24) 



which is related to the absorption coefficient a of a 
Doppler-broadened medium by 



a = T] (im 



P3l(^) 



(25) 



The absorption gives rise to the attenuation of the prop- 
agating light pulse, that is described by the Beer's law, 
namely. 



R21 {t,T) 
R21 (0,T) 



-2al 



-2aVat 



(26) 



where / is the propagation distance and Vg = ^l/^r] is 
the group velocity. It is straightforward to obtain hz = 
4:{kvs/^c)'^ T with eq.(26). In contrast to the Gaussian 
decay in eq.(22), which depends on temperature only, the 



decay of ground-state coherence in SL depends on the 
temperature, coupling field and the spontaneous decay 
rate of level |3). The numerical simulations of R21 (t^T) 
and the analytic predictions are plotted for various tem- 
peratures in Fig. 6, where good agreements are demon- 
strated. It is not unexpected that our numerical sim- 
ulations closely agree with those obtained by averaging 
the solutions of optical-Bloch equations subjected to the 
Maxwell-Boltzmann distribution at a particular temper- 
ature, since in our numerical calculations so far, the effect 

of recoil is negligible, i.e., — Vg ^ h{kp — kcY /2m 



we have A„ = kv and 



-kv in the above equation. 



[see the definition of below eq.( 20 )] . However, we note 
that the above criterion is no longer valid if the mass of 
atom is made small and the coupling and probe beams 
are counter-propagating. The resultant dephasing can 
significantly reduce the output level of the probe pulse 
and our scheme appears applicable to this kind of prob- 
lems. 

Finally, it should be noted that, simply by adjust- 
ing the temperature (or equivalently, the velocity dis- 
tribution), we can simulate the decaying behavior of the 
ground-state coherence, provided that the explicit time- 
dependence of the coupling field is given. This can not be 
achieved via solving the optical-Bloch equation by impos- 
ing a phenomenological decay rate 7 on the metastable 
ground state |2) (which is accurate only at Vtc ^ iii 
SL and not valid at all in LS), and thus features the ma- 
jor decoherence between our formalism and the optical- 
Bloch equation. For this reason, it is expected that the 
current scheme can be applied to investigate the dynam- 
ics of stationary light pulse (SLP) which is basically con- 
sisted of four SL processes — two co-propagating and two 
counter-propagating. Owing to the inherent complexity, 
specifying the high-order 7's in the SLP turns out to be 
tricky when solving the usual optical Bloch equations [25] 
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. Under the circumstances, without doubt, our scheme 
appears to be an easier and more natural way to approach 
the problem. 
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FIG. 5: (Color online) The triangles, squares and circles are 
the numerical results of the normalized i?2i (t, T) during the 
storage at T = 100, 300 and 500 /xK, respectively; the dashed- 
dotted, dotted and dashed lines are the predictions by eq.(22 ). 



IV. CONCLUDING REMARKS 



We have presented a numerical scheme to study the 
dynamics of SL and LS in an EIT medium at finite tem- 
peratures. Based on the gauge invariance of Schrodinger 
equation under Galilean transformation, we derive a set 
of coupled equations for a boosted closed 3-level EIT sys- 
tems. The loss of ground-state coherence at finite tem- 
peratures is then treated as a consequence of superpo- 
sition of density matrices representing the EIT systems 
moving at different velocities. Unlike other theoretical 
treatments in which atoms are assumed immobile, our 
scheme takes both atom's external and internal degrees of 
freedom into full account. The feasibility of this scheme 
is shown by comparing the numerical results to the ex- 
perimental data for both SL and LS. Last but not least, 
this scheme also enables us to study the dynamical prop- 
erties of a Doppler-broadened EIT medium in the non- 
perburbative regime of probe and coupling pulses with 
comparable intensities, in which new effects are expected 
to arise. 
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FIG. 6: (Color online) The normalized R21 {t,T) with Qc — 
0.75r, OD = 400; the triangles, squares and circles are the 
numerical results at T = 100, 300 and 500/iK, respectively; 
the dashed-dotted, dotted and dashed lines are the predictions 
of eq.(l26l. 
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